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In the last years there has been an increasing interest in the study of the current corre- 
lators in QCD at finite temperatures. The hope is that investigating the same correlators, 
both at high temperatures, where the state of quark-gluon plasma is expected, and at low 
temperatures, where the hadronic phase persists, a clear signal for a phase transition could 
be found. For a review of calculations of correlators performed by various analytical methods 
and in the lattice simulations, see, e.g., Ref. 

The study of temperature dependence of current correlators is interesting in many as- 
pects. At small distances the correlators are expressed through the operator product expan- 
sion (OPE) in terms of matrix elements of the operators of low dimension. In this way, the 
temperature dependence of these matrix elements manifests itself in the temperature depen- 
dence of the correlators and visa verse. At T = using dispersion relations the correlators 
can be expressed in terms of contributions of hadronic states. Then, using some theoretical 
tools (differentiation, Borel transformation, etc) it is possible to enhance the contribution of 
lowest hadronic states. Therefore, knowledge of T-dependence of hadronic correlators can 
give us an information about the T-dependence of masses of the lowest hadronic states. Since 
in the approach of QCD sum rules, these masses are determined by the matrix elements of 
operators in OPE, the T-dependences of both are interrelated. 

As argued in Ref. f^, in carrying out this program certain wrong steps had been taken, 
and a misunderstanding existed in the literature. Thus, some authors maintained that the 
thermal shift of hadron masses occurs already in the lowest order in temperature, O(T^). 
However, the general statement is that hadron masses do not move in this order p|, This 
statement is very general: it is based only on PCAC and current algebra. (The result that 
the nucleon pole does not move in order was obtained in the chiral perturbation theory 
in Ref. Q and by considering a current correlator in Ref. 0). The only interesting physical 
phenomenon that occurs in this order is the parity mixing, i.e. the appearance of states with 
opposite parity in a given channel and, in some cases, also an isospin mixing. The latter arises 
in baryon current correlators where, for example, in the current with the quantum numbers 
of A there appears a S pole, and in the nucleon channel there appear poles corresponding 
to baryon resonances with = and T = |, |. For the case of vector and axial currents 
in the chiral limit, this mixing is given by 



where C^^^\q^T) are the correlators of V and A currents at finite temperature, CY^I;^\q,Q) 
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are the same correlators at T = 0, 



and -Ftt = 93 MeV is the pion decay constant. 



(2) 



At T 7^ both and Cf^^, have transverse and longitudinal parts 



Cl!;^\q.T) = i-g,^q' + q,q.)Cr^^\q\T) + q,q.C^^^\q' ,T) (3) 



At T = CY {q^, 0) = 0, but C(^{q^^ 0) is non-zero and in the chiral limit (when light quarks 
and pions are massless) is given by the one-pion contribution 



Cf{q\Q) = ^ (4) 



According to Eq. at T 7^ the longitudinal part (the pion pole) appears also in the 
vector channel. 

If C^'^(g, 0) are represented through dispersion relations by contributions of the physical 
states in the V and A channels (p, ai), then according to Eq. (|l|) the poles that are in the 
r.h.s. of Eq. (|l|), i.e. at T = 0, appear at the same positions in the l.h.s. Therefore, in order 

the poles corresponding to p and ax do not move p[ . An important consequence of Eq. (||) 
is that at T 7^ in the vector (transverse) channel apart from the poles corresponding to 
vector particles, there arise poles corresponding to axial particles and vice versa, i.e. a sort 
of parity mixing phenomenon occurs. The manifestation of this phenomenon is in complete 
accord with the general considerations presented in Ref. 0: the appearance of the p (ai) 
pole in the axial (vector) channel corresponds to singularities in the s-channel. In the same 
way a pion pole appears in the longitudinal part of the vector channel. 

In this paper we show that in the next order, O(T^), such a simple picture, where 
the current correlator at finite temperature is represented by the superposition of T = 
correlators, does not take place. Interpreted in terms of temperature-dependent poles, it 
would mean that masses are shifted in this order. In what follows we consider only the 
transverse part of the correlator because the T-dependence of the pion mass and the decay 
constant was thoroughly investigated earlier 0, and we can say nothing new here. 

The thermal correlation functions in Eq. (|I]) are defined as 





|2/c?W^-T{j^(a;),j,(0)}e-^/'^| 


\n) 






\n) 



C,M,T) = {t d'xe^^^T{Ux),m})T = ^^^^^^ / Y-HT \ (5) 
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where H is the QCD Hamiltonian, the sum is over all states of the spectrum, and j^{x) is 
either a vector or an axial current, i^{x) = (1/2)('U7^(75)m — d'^^{;^^)d). It is assumed that 
is space-like, = —q^ > 0, and is much larger than a characteristic hadronic scale, 
^ R~^, where Rc is the confinement radius, ~ 0.5 GeV. 

We consider the case of temperatures T below the phase transition temperature T^. In 
principle, the summation over n in Eq. (|^) can be performed over any complete set of states 
\n) in the Hilbert space. It is clear however, that at T < Tc the suitable set of states is 
the set of hadronic states, but not the quark-gluon basis. Indeed, in this case the original 
particles forming the heat bath, which is probed by external currents, are hadrons. The 
summation over the quark-gluon basis of states would require to take into account the full 
range of their interaction. In connection with consideration of current correlators at finite 
T, this point was first explicitly made in Ref. [Q]. In the early papers [0 devoted to the 
extension of QCD sum rules to finite temperatures the summation over |n) at low T was 
performed in the quark-gluon basis without account of confinement. In a recent paper 
QCD sum rules at finite T were rewritten using the pion basis. 

At T well below the phase transition temperature an expansion in T can be performed. 
The main contribution comes from the pion states, \n) = |7r), |27r), ... In this paper we 
restrict ourselves to the chiral limit, when u, d quarks and pions are massless. The corrections 
to the chiral limit will be considered in a separate publication. 

In the chiral limit there are two parameters in the low T expansion. One parameter 
appears when the pion momenta, p ~ T, in the matrix elements in Eq. can be neglected. 
Then, the matrix elements in Eq. can be calculated using PCAC and current algebra. 
The powers of arise due to phase space integration with the Bose factor. In this case the 
expansion parameter is T'^/F^: the one-pion contribution is proportional to T'^/F^, the two- 
pion contribution is of order T'^/F^, etc 0. In the order there are only terms of this type. 
It is clear that in any order (T^)^ the terms of this type are expressed through the vector and 
axial vacuum correlators C^j'^\q'^,0) and, as a consequence, do not result in thermal shifts 
of hadron masses. The other expansion parameter appears if non-vanishing pion momenta 
in the matrix elements in Eq. (^) are taken into account. Since the characteristic distances 
in Eq. (1) are of order ~ l/Q^? the expansion parameter in this case is T^/Q^. It is also 
worth to mention that the contributions of massive hadronic states \n) are exponentially 
suppressed as exp(— m„/T). 

Let us start with the calculation of terms of the first type. The matrix elements 
over two-pion states in the limit p give a contribution. But this is not just a 
second iteration of the procedure used to obtain the one-pion matrix element. It is also 
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necessary to take into account the interaction between the pions in the initial and the finite 
states. This can be illustrated by the example of terms in the T-dependence of the quark 
condensate 0: 

mr = {qq)o (l " " ^^') (6) 

The e and terms here come from one- and two-pion matrix elements, respectively. However, 
not accounting for the initial (finite) state interaction of pions in the two-pion state would 
give 9e^/32 instead. The interaction amplitude for tt^tt'^ —>■ tt^tt'^ to zeroth order in pion 
momenta is given by 

-T^i^SacSbd - SabScd - SadSbc) ■ (7) 

The in the numerator of the above amplitude cancels against the pion mass in the pion 
propagator at p — > and contributes — 3e^/8 to the T-dependence of (qq), making the correct 
total of -3eV32. 

Similarly, it is easy to show that in the case of vector and axial correlators the two-pion 
matrix elements with the account of initial (finite) state interaction of pions result in the 
following expressions for the transverse part of V and A correlators: 

C^}^\q,T) = {-g,.q' + q,q.) {Cy^A){q') + {e - ^6^) (C^(v)(g^) - Cv^A){q'))) • (8) 

As was expected, the T^/F^ terms above are again expressed through correlators at T = 0, 
and they do not result in a thermal mass shift. Also, these terms are Lorentz invariant, since 
only the tensor —g^uq^ + q^J.qv appears here. 

Consider now the terms of the second type, that arise from non-zero pion momenta. It is 
well known from the description of deep inelastic electron-hadron scattering Q that a matrix 
element of the product of two vector currents may be represented using two tensor structures 

T;Av.q)=^ J dV^iT{p)\T{j,{x),UO)}Kp)) 

= [-9,. + '-f) T,(., g^) + [p, - [p, - T.(., g^) , (9) 

^We are now considering the case of a vector correlator. However, the axial correlator in the chiral limit 
has the same tensor structure, and the results obtained in the vector case may be directly applied to the 
axial one. 
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where v = pq, Ti is dimensionless, and T2 has dimension mass~^. The corresponding 
contribution to the thermal correlation function is obtained by integrating the above equation 
over the pion phase space with the Bose factor. And it is the second term in the r.h.s. of 
Eq. (j^) that, after this integration, provides the expected Lorentz non-invariant contribution 
to the thermal correlator. In terms of the OPE, the function T2 is contributed only by 
averages of Lorentz non-scalar operators, while Ti receives contributions from both Lorentz 
scalar and non-scalar operators]^. 



It is well known [|10| that the function T2 satisfies a dispersion relation without subtrac- 
tions: 

Having in mind the subsequent integration over we are interested in T2(0, g^). Then, using 
the relation between the imaginary part of T2 and the structure function 

ImT2(z.,g2) = -F2(x,g2), (11) 
where x = (5^/2z/ is the standard deep inelastic scaling variable, we get 

T2{0,q') = -^ J\2ix,q^)dx. (12) 
Similarly, the function Ti(z/, g^) satisfies a dispersion relation with one subtraction, 

T,{u,q')-T^{0,q') = — U'^ 2, • 13 

Again, using the relation between ImTi and the structure function Fi 

Imri(x,g2) =27rFi(x,g2), (14) 

we obtain 

Ti(i/, q^) - ri(0, q^) = ^ £ 2xF,ix, q')dx . (15) 

The subtraction constant Ti(0, g^) corresponds to zero momenta of the initial and final pions 
and was already accounted for in the terms proportional to T'^/F^. It is worth to mention 
that at this stage the assumption » was not used. In the derivation of Eqs. (|12D 
and (|15D it was assumed that u << Q^, which is equivalent in the chiral limit to Q > T. 
The same assumption was sufficient for the derivation of Eq. (H). Therefore, Eqs. (P), (|12D 
and (|T5|) are correct also at moderate ~ 1 GeV^. 
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At higher Q^, in the scahng region, the integrals in Eqs. (|T2D and (|T5D are equal, since in 
this region 2xFi{x,Q'^) = F2{x,Q'^). The integral 



M2 = F2{x)dx (16) 

JO 

is the second moment of the structure function and in the parton model has the meaning of 
the fraction of the pion momentum carried by quarks. In our normalization of currents 
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1 1 /■! 1 /■! 

-■-/ X dx y][q^a{x) + q^a{x)] = - x dx[v{x) + 2s{x)], (17) 
4 3 JO 2 JO 

where v{x) and s{x) are the distributions of valence and sea quarks in the pion (see, e.g. 
p!0|). The factor 1/4 in Eq. ([T7|) comes from the definition of the currents in Eq. (^, and 
the factor 1/3 arises due to averaging over vr"'", vr", 7r° in the heat bath. 



The distributions v{x) and s{x) were parametrized in Ref. to fit the experimental 
data on the Drell-Yan process vr + — > + X and on the direct photon production 
77 + — > 7 + X, and it was found that M2 ~ 0.12, which is somewhat lower than the 



result 0.15 ± 0.02 of Ref. |]lT|, where this quantity was obtained from QCD sum rules for a 
correlation function in external symmetric tensor field, and close to 0.11 ±0.03 as estimated 
in Ref. ||^. These numbers correspond to the normalization point /i = 1 GeV. 

Now, to obtain the corresponding contribution to the thermal correlator, one has just to 
do the Bose-weighted integrals over the pion momentum and sum over the three pions 

'^'^ ^ t:m<i) (18) 



{2iif2p exp(|Mp|/T) - 1 
[u is the 4- velocity of the heat bath), which gives together with Eq. (§) 

C,u{q, T) = (-g^^q^ + q^q,) C,{T, q') + «C2(r, q') , (19) 
where u^^ = u^- {uq)q^/q^, 



C^{T,q') = Cy{0, q') +(e- \e^) (C^(g^) - Cy{q')) + 1 c ^ + 



and 



q^ 

C2{T,q^) = -c^, (20) 



c^^. (21) 
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The low temperature expansion of Ci contains powers of both T'^/F^ and T'^/Q'^, but powers 
of T'^/F^ are absent in C2. Notice, that while all three pion charge states contribute to T^/Q'^ 
terms, only two of them contribute to T'^/F^ terms. The same formulae hold for the axial 
correlator, with the obvious change V A. 

The above formulae for the thermal correlator may be considered also from the viewpoint 
of the OPE for the correlation function(see, Refs. H, ^). The OPE itself of course carries no 
information hey about the state over which the matrix element of the operators is considered, 
or about the heat bath, in case of finite T. It contains operators of arbitrary Lorentz spin 
s and twist t. When vacuum correlators are considered, only scalar, s = 0, operators 
contribute, while non-scalar, s ^ 0, operators drop out in the averaging. When averaging 
over a hadron state or over a heat bath, the s 7^ operators do contribute, since there is an 
additional vector in the problem, the target momentum p or the 4- velocity of the heat bath 
u. The matrix elements of s = operators over pions or the heat bath may be estimated (if 
the pion momenta can be neglected) using PCAC, which relates them to vacuum averages. 
The s = operators do not give (T^/Q^)" terms and therefore do not result in thermal 
mass shifts. The general expression for the hadron matrix element of an s = ra, n > 2 
operator is (p|0^i,/i2,.--Mnb) = '^PAtiPM2---PMn (i^ ^^e chiral limit there are no trace terms 
~ Ofitfij)- These matrix elements cannot be reduced by PCAC to vacuum averages and are 
new non-perturbative parameters. Their Bose-weighted integrals over the pion momenta are 
T-dependent s 7^ condensates which are suppressed as compared to the T-dependent 
parts of s = condensates. 

It is clear that in terms of the OPE the function Ci is contributed both by s = and 
s 7^ condensates, as is Ti. However, C2 and T2 are related only to the s 7^ condensates. 

In the chiral limit a difference in the s = operators for the vector and axial correlators 
appears on the level of 4-quark operators 0. A good consistency check of the calculation of 
correlators in the pion gas approximation is on whether the T dependences of the s = 
4-quark condensates in C]"^^^ match the V — A mixing in the first of Eqs. (^). This indeed 
turned out to be the demonstrated (to order T^) in Ref. |T^, and also in Ref. [|T3| for 



baryonic currents. The correlation in T-dependences of the s = condensates in opposite 
parity channels is not accidental and is related to the scattering of thermal pions on the 
currents. Notice also, that since 

C!'(T) - C^{T) = (1 - 2e + e^) (^^(0) - C^{0)) (22) 



^The gluon condensate gets its T dependence only in order anyway. 
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this correlation exactly satisfies Weinberg sum rules generalized to finite T. 



Among the non-singlet condensates, the leading contribution to C2 at low T comes from 
the lowest spin, s = 2, which corresponds to the behavior. In the leading twist there are 
two s = 2 operators which are related to the energy- momentum tensors of quarks, 61^, and 
gluons, e^^, 



1 

= ^{QlxD^q + q-faDxq), q = u , d 



Explicit expressions for the contribution of these operators to Ti and T2 can be obtained 



from the general formulae of the theory of deep inelastic scattering (see, e.g. |17]). We present 
here the result for the case, when the longitudinal structure function Fi = 2xFi{x) — F2{x) 
is neglected and only QCD corrections proportional to lii(Q^//i^) are retained. It can 
be shown that since all pion charge states are equally populated in the heat bath, only 
flavour-singlet operators contribute to the structure functions. Then 



Uu,q') = ^T2{u,q') 
2x 



q\qo 

q' 



1 



(24) 



where the averaging over the three pion charge states must be performed. It is easy to see 
that the contribution of 9'^^ to Eq. (^) is small at ~ 1 GeV^. 

The pion matrix element of the total energy- momentum tensor is just a normalization 
constant, 



{{n{p)\n{p')) = {2nf2E5'^^\p - p')). On the other hand, we have 



(25) 



{Ap)K + OiMp)) = 8M2p,p. . (26) 

So, if we define also a constant b as (vr|6'^j^|7r) = bp^p^, then 8M2 + 6 = 2. The constant b 
enters the matrix element (7r|E^ + B^lvr) and also (E^ + B^)^. It was determined in Ref. 
from a duality type of consideration: b = 1.16 ±0.14 at fi = 1 GeV. This is in accord with the 
estimates for M2 obtained in Refs. M, |ll| and with the statement that gluons carry about 



50of nucleons. 
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Now, we would like to discuss a possibility to interpret the T'^/Q'^ corrections to the 
correlators in Eqs. (^) in terms of particle thermal mass shifts. It is convenient to use the 
standard representation [|^ of the vector correlator in a medium in terms of two invariant 
functions Ci and Ct-, which in the rest frame of the heat bath {u = (1,0)) are defined as 

Cl = {^,-'-f)Cl+'-f<llCf, (27) 

or in terms of Ci and C2 

Ct = q'C, 

Ci = C, + ^C2. (28) 
Q 

At T = these two functions are not independent, Ci = Ct/q"^ = Cy. They are also related 
at T ^ 0, if 

q = 0, go ^ 0, = -gg 

cYiT) - cYiT = 0) = cr - cr' = {cj - cj-') 

--{^-V)^CX-Ct)^cI^,. (29) 

where c was defined in Eq. (pTf). 
In another special case 

go = 0, q ^ 0, = 

cj-cr' = -{^-r')(Cy-c^) + ^^2^ (30) 

Till now we considered the correlators at negative g^ = —Q^. In order to interpret the 
results in terms of the particle mass shifts, the amplitudes at negative g^ must be represented 
using dispersion relations through the contributions of physical states, defined at positive 
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and s. Unlike the case of T = 0, where the correlators are functions of one variable, 
g^, at T 7^ they are functions of two variables, go and q^. In this case the only way to 
represent the amphtude at negative through the contributions of physical states is to use 
the dispersion relation in go at fixed q^. (In the opposite case, when go is fixed and q^ is 
variable, the amplitudes would have non-physical singularities). 

So, let us consider the case q = 0, go 7^ 0, = —q^. In this case there is only one 
structure function, C/(go) = Ci(go). We choose the standard model for the spectral densities 
as a sum of the lowest resonances and continuum. The dispersion relations over go are 
contributed by the physical states in the g^- and s-channels. Therefore, for the structure 
function C^iQ'^.T) we have 

Cr(Q^T)^J^.J^.'-r'-^i^. (31) 



where the first and the second terms in the r.h.s. of Eq. (^) correspond to the contributions 
of p and ai mesons to the vector current correlator, Ap y(T) and A^^ y(T) are the correspond- 
ing T-dependent coupling constants, and Xai,v ~ T'^- (The subtraction constant is omitted). 
A similar equation holds for the axial current correlator Cj^((5^, T). From Eq. (^) it is easy 
to see that the terms of order T'^/Q^ vanish in the difference CY{Q^,T) — C^{Q^,T): 



e 
2 



-e l 



(32) 



If we put 



mJ(T) =ml + 6ml, m\ (T) = m\ + bm^ 



a\ 1 



then from Eqs. (^) and (|32|) it follows that 



(33) 



(34) 



and 



\2 



(T) 



X2 



(T) 



l-e 1 



e l 



(35) 
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For the sum CY{Q^ T) + C^^(g^ T) from Eq. (H) we have 

CX{Q\T) + Ct{Q\T)- 

It is clear from the comparison of Eqs. (^) and (^) that with our model of hadronic 
spectrum the continuum cannot contribute to the l.h.s. of Eq. (P^), since the imaginary part 
of the r.h.s. vanishes at < 0. Then from Eqs. 



cnQ\o) + c^{Q',o) 



n2 



(36) 



) and (Bq) we have 



Xl5ml + Xl6ml=-cT\ (37) 

(Taking into account only terms ~ we put X^iT) ^ A^(0) = and Xl^{T) ^ Xl^{0) = 
XlJ. From Eqs. (|3|) and (0) we have 

'^"^?=-c^' '^"^ai = -c^. (38) 

The residues in Eq. (P5| ) are the standard couplings of p and ai mesons with the vector and 
axial currents, A^ = m'j^/gp, A^^ = 'm'^_^/gl_^. Numerically they are rather close, A^ ~ A^^ ~ 
0.02 i. 

We see that both the p and ai masses start decreasing with T, and the mass shifts appear 
in order and, in terms of OPE, are due to Lorentz non-scalar condensates as emphasized 
in refs.0, |14| 0. The corrections proportional to powers of T'^/F^ affect only the residues 
of the currents. This fact can be easily understood. Indeed, in the OPE for the correlators 
taking into account finite temperatures to order would result only in the change of the 
same Lorentz scalar condensates which appear in OPE at T = 0. Then it is clear from 
the representation through dispersion relations that any such change can be described by 
modifications of the residues without affecting the position of poles. 

The scenario when both the vector and axial masses decrease with T is allowed by 

Weinberg sum rules at T 7^ 0|T^. Numerically the mass shifts are rather small. Even 

at T = 150 200 MeV (usually accepted values for the phase transition temperature) 

(5m^ ^ 6m?^_^ ~ —(0.01 — 0.02) GeV^. At the same time at these temperatures the change of 

residues according to Eq. ( ^5| ) is very essential. (For this reason Eqs. (|38|) are not completely 

reliable at T > 100 MeV, since we put X^{T) ^ A2(0)). 

•^In Ref. it was claimed that the p and ai masses move in order T^. This conclusion comes from an 
incorrect choice of the spectral density in the sum rules: the mixing of vector and axial channels was not 
taken into account. For this reason the results obtained in this paper for the T-dependence of the masses 
are not reliable, though we agree with the analysis of T-dependent condensates carried out there. 
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Let us summarize our main results. The corrections of order T"^ to the correlators of 
vector and axial currents were calculated in QCD in a model independent way in the chiral 
approximation when u, d quarks and pions are massless. The results are expressed in terms 
of the second moment of the pion structure function in deep inelastic lepton-pion scattering 
that is equal to the matrix element of the quark energy-momentum tensor over the pion state, 
or to the fraction of the pion momentum carried by quarks in the parton model. Interpreted 
in terms of physical mesons the calculated corrections correspond to negative mass shifts 
of the p and ai mesons proportional to T^. (As was shown earlier P, the mass shifts 
are absent in order T^). These mass shift originate from Lorentz non-scalar condensates in 
OPE. Numerically they are rather small at T < 100 MeV, where our approach is correct. 
The corrections arising from a finite pion mass were not touched in this paper. We plan to 
consider them in a future work. 
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ported in part by the International Science Foundation under Grant No. M9II000. The 
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